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ABSTRACT 


Calculations  have  been  made  of  "flame"  (or  mixing)  dimensions 
associated  with  afterburning  in  rocket  exhaust  at  high  altitudes  for 
a  hypothetical  rocket  and  trajectory.  From  these  and  a  given  exhaust 
composition,  one  can  estimate  concentrations  and  times  available  for 
reactions  of  interest.  Above  some  altitudes  if  observed  flame  lengths 
should  turn  out  to  be  much  greater,  it  would  indicate  that  reaction 
rates  rather  than  mixing  rates  were  controlling. 

Investigation  of  negative  ion  formation  reactions  indicates  that 
they  may  not  be  a  significant  factor  in  removal  of  electrons,  since 
dissociative  recombination  of  polyatomic  ions  is  much  faster.  However, 
if  the  production  of  electrons  were  mainly  due  to  ionization  of  alkali 
metals  (which  have  very  low  recombination  coefficients),  the  above 
reactions  could  be  an  important  electron-removal  mechanism. 

Calculations  have  been  made  of  electron  densities  resulting  from 
thermal  ionization  of  potassium  in  a  rocket  engine  at  various  altitudes 
to  provide  estimates  for  comparison  and  evaluation  of  the  probable 
significance  of  other  reactions.  The  estimates  indicate  that  thermal 
ionization  of  potassium  with  typical  solid  propellants  can  be  a  very 
important  source  of  electrons. 

Preliminary  investigation  of  photoionization  in  rocket  exhaust 
has  provided  some  values  of  radiation  flux  density  versus  altitude  tor 
various  wavelengths. 


Reviewed  and  approved: 
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INTRODUCTION 


Free  electrons  In  a  rocket  exhaust  stream  cause  Interference  In 
communications  with  the  vehicle  and  Increase  Its  vulnerability  to  de¬ 
tection.  An  important  source  of  free  electrons  is  the  chemionlzation 
produced  in  the  chemical  reactions  of  afterburning  as  unburned  fuel  in 
the  exhaust  mixes  with  the  ambient  atmosphere. 

This  study,  which  is  a  continuation  of  the  work  performed  under 
Contract  AF  04(647)-751,  surveys  existing  knowledge  of  afterburning 
phenomena  in  rocket  exhausts  related  to  the  production  of  ions.  The 
objective  is  to  provide  a  summary  of  existing  information  and  to  suggest 
a  program  of  theoretical  and  experimental  research  to  fill  gaps  in 
present  knowledge.  The  over-all  aim  is  to  predict  or  control  such 
ionization. 

Efforts  will  be  concentrated  on  two  general  areas:  (l)  fluid 
dynamics  of  exhaust  jets  in  the  atmosphere  at  various  altitudes,  and 
(2)  ionization  reactions  which  may  occur  during  afterburning.  The  two 
areas  are  interrelated  and  cannot  be  arbitrarily  separated.  However, 
because  of  the  complexity  of  each  field  it  is  desirable  to  treat  them 
independently  insofar  as  possible. 

During  this  report  period,  flame  length  calculations  presented  in 
the  First  Quarterly  Report*  have  been  treated  in  more  detail,  and  studies 
have  been  made  of  (a)  the  possible  significance  of  negative  ion  forma¬ 
tion  as  electron  sinks,  (o)  the  importance  of  alkali  metal  contaminants 
in  the  propellant,  especially  in  relation  to  other  electron  sources,  and 
(c)  the  possible  role  of  photoionization  in  rocket  exhaust. 

Contributions  to  this  report  have  been  made  by  several  people.  The 
following  sections  are  primarily  the  work  of  the  following  individuals: 

Negative  Ion  Formation  —  Robert  C.  Whitten 

Thermal  Ionization  of  Alkali  Metals  —  Carole  R.  Gatz 

Photoionization  —  I.  G.  Poppoff 
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DISCUSSION 


A  Flame  Length 

The  fluid  dynamics  of  an  exhaust  jet  is  extremely  complicated. 
However,  at  high  altitudes  certain  simplifications  are  admissible.  For 
example,  shockwave  structures  within  the  jet  can  be  neglected  at  high 
altitudes  and  high  velocities . 2 > 3  Further,  at  great  altitudes3  (exceed¬ 
ing  60  km)  the  expansion  of  the  jet  to  roughly  a  cylindrical  shape  is 
very  rapid,  so  that  the  boundary  of  the  jet  can  be  approximated  by  a 
blunt  body  with  a  cylindrical  afterbody.  Thus,  as  an  approximation 
the  jet  is  assumed  to  expand  to  a  cylindrical  shape  at  a  distance  b 
from  the  nozzle  exit,  and  mixing  is  assumed  to  begin  at  x  =  0  (see 
Fig.  1).  A  shock  will  be  caused  by  the  jet  "blunt  body,"  and  as  the 
diameter  of  the  expanded  jet  increases  (at  higher  altitudes  due  to 
lower  ambient  pressure),  the  shock  stand-off  distance  may  actually  ex¬ 
ceed  the  length  of  the  rocket,  and  thus  the  bow  shock  from  the  expanded 
exhaust  jet  can  occur  ahead  of  the  rocket  Itself.  Shock  stand-off 


FIG.  1  SCHEMATIC  REPRESENTATION  OF  ROCKET  EXHAUST  JET 
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distance  A  for  a  sphere  can  be  calculated4  from 

A/'r0  =  0.143  exp  (3.24/Mj) 

where  R0  =  radius  of  the  sphere,  and  =  free-stream  Mach  number.  For 
our  region  of  interest,  A  will  be  roughly  0.15  times  the  radius  of  the 
expanded  Jet.  For  our  simple  model,  we  will  neglect  the  effects  of  the 
rocket  body,  since  the  expanded  jet  will  be  large  in  comparison  with  it. 

The  pressure  at  the  stagnation  point8  of  a  hemisphere  is  assumed 
to  be  given  by 


p  ,,/p  =  7  M*  +  1 

*s  t  CC  IX)  00 

where  p  =  free-stream  pressure  and  M  =  free-stream  Mach  number.  The 

00  00 

pressure  decays  around  the  hemisphere  toward  the  cylinder;  at  the 
transition,  corresponding  to  x  =  0  (Fig.  l),  the  pressure  is  assumed®*7 
to  be  0.1  pgt .  Decay  to  ambient  pressure  may  be  assumed  to  occur  at 
25-100  diameters  behind  a  sphere,®  or  at  about  150  diameters  on  a 
cylindrical  afterbody . ® > 8  In  the  model  for  the  first  calculations  of 
flame  length1  it  was  assumed  that  ambient  pressure  was  attained  at 
x  =  0  (an  unspecified  distance  b  from  the  front  of  the  expanding  exhaust 
plume)  and  mixing  was  assumed  to  begin  at  this  point.  However,  if 
x  =  0  is  at  about  100  expanded  jet  diameters  downstream,  one  should  assume 
that  the  disturbance  which  initiates  the  growth  of  turbulent  mixing  occurs 
at  an  earlier  point. 

For  our  present  model,  the  distance  b  to  where  x  =  0  is  assumed  to 
be  roughly  one  jet  radius,  i.e.,  the  blunt  portion  of  the  jet  is  assumed 
to  be  roughly  hemispherical,  the  pressure  at  that  point  is  taken  as 
0.1  Pgt i  and  mixing  is  assumed  to  begin  there.  In  order  to  obtain  some 
dimensions  for  a  jet  flow  field  at  high  altitudes,  calculations  by 
W.T.  Lord10  for  a  particular  rocket  and  specific  trajectory  have  been 
used.  The  semi  -dive  rgence  angle  of  the  nozzle  is  15°,  the  ratio  of  the 
specific  heats  of  the  jet  gas  is  1.2,  the  velocity  of  the  jet  gas  at 
the  exit  plane  is  3000  meters/ sec,  and  the  pressure  is  1/3  atmosphere. 
Diameter,  2a,  of  the  expanded  jet  is  presented  for  different  altitudes 
(Table  I,  and  Fig.  2).  The  rocket  velocity  varies  linearly  with  altitude 
and  reaches  a  value  of  7300  meters/ sec  at  150  km. 
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|  FIG.  2  EXPANDED  EXHAUST  JET  DIMENSIONS  (Reference  10) 

Expanded  Jet  velocity  (average)  is  obtained  from  the  relation  for 
conservation  of  mass  flux. 

» 

°j  “j  *j  ‘  " 

where  =  density,  =  velocity,  and  =  cross  sectional  area,  each 

referring  to  the  expanded  Jet,  and  M  is  the  mass  flowrate  at  the  nozzle. 

Jet  temperature  was  calculated  assuming  isentropic  expansion.  Ambient 

temperatures,  pressures,  densities,  etc.  were  taken  from  the  Handbook 

of  Geophysics.11  External  air  temperature,  T  ,  refers  to  the  air  within 

e 

a  few  jet  plume  radii,  which  has  been  shock-heated  in  passing  through 
the  bow  shock,  and  has  been  obtained  from  calculations  by  Hochstim1*  for 
equilibrium  air.  External  air  velocity,  ug  (adjacent  to  the  jet)  is 
assumed  to  be  0.5  -  0.75  times  free-stream  velocity,  u  .  For  a  sphere, 
Feldman8  quotes  a  value  of  0.75  u  at  the  downstream  distance  where 

I  00 

;  expansion  to  free-stream  pressure  has  taken  place.  For  purposes  of  our 

calculation,  we  have  taken  u  =  0.75  u  . 

e  oo 

In  order  to  make  some  specific  calculations  of  flame  lengths  and 
diameters  applicable  to  a  rocket  rising  through  the  atmosphere,  a  hypo- 
I  thetical  rocket  and  trajectory  have  been  defined. 10> 1  *  Exhaust  compo- 

|  sition1 *  is  given  in  Table  II. 
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Table  II 


ROCKET  EXHAUST  COMPOSITION 


Constituent 

CO 

Ha0  H2 

C02 

OH  H  02  0 

Mass  %  of  Total 

46 

27  1 

23 

2  0.1  0.6  0.3 

Mass  ratio  air/exhaust  required  to  complete  combustion  to 
C0a  and  Ha0  —  1.55 

Average  molecular  weight  of  exhaust  —  22.1 
Exhaust  exit  temperature  —  1800 °K 


A  calculation  of  the  flame  shape  in  a  jet  formed  by  a  uniform  cylin¬ 
drical  stream  of  hydrogen  surrounded  by  an  ambient  parallel  high-speed 
flow  of  air  at  a  pressure  equal  to  that  of  the  jet  (see  Fig.  3)  has  been 
published.14  This  calculation  is  based  on  a  theoretical  analysis  involv¬ 
ing  a  hypothetical  model15  for  the  eddy  viscosity  relating  it  to  the  eddy 
viscosity  for  incompressible  flow.  The  published  plot  of  results  of  the 
calculation1 4  is  presented  in  Fig.  4.  Here  rf  is  the  radius  of  the 
flame  surface  at  the  axial  distance  x  from  the  plane  where  mixing  starts 

(l.e.,  where  x  =  0);  a  is  the  radius  of  the  jet  at  x  =  0;  and  u./u  is 

j  e 


FIG.  3  SCHEMATIC  REPRESENTATION  OF  FLOW  (bo.ed  on  Roforonco  14) 
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x/2a 


HA- 4134-5 


FIG.  4  SHAPE  OF  FLAME  SHEET  (Reference  14,  Fig-  10) 


the  ratio  of  jet  velocity  to  external  flow  velocity,  relative  to  the 
nozzle.  In  Fig.  4  it  is  assumed  that  the  chemical  reaction  rate  is 
infinite  and  that  the  turbulent  flow  and  mixing  process  control  the 
flame  surface  location. 

Flame  lengths  and  maximum  diameters  can  be  obtained  from  Fig.  5 

(derived  from  Fig.  4)  as  a  function  of  u  /u  .  Flame  length  L  (or  L/2a 

J  ® 

in  dimensionless  form)  is  defined  as  the  value  of  x  (or  x/2a)  when  rf/a 
is  zero.  However,  these  values  refer  to  the  following  special  conditions 
a  jet  of  hydrogen  gas  at  1300°K,  external  flow  of  air  at  1600°K,  both 
at  1-atmosphere  pressure,  and  ambient  air  temperature  (far  from  the 
exhaust  flow  field)  of  282°K.  Application  of  Libby’s  method  to  other 
gases  and  conditions  (different  molecular  weights,  concentrations,  tem¬ 
peratures,  etc.)  requires  very  lengthy  calculations,  unless  one  can  use 
some  simple  theory  of  mixing  and  assume  some  relation  between  it  and 
Libby's  theory.  This  has  been  done  to  provide  a  preliminary  test  of 
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FIG.  5  FLAME  DIMENSIONS  (derived  from  Fig.  4) 


Libby's  method.  Use  has  been  made  of  some  work  by  Ricou  and  Spalding16 
on  determination  of  entrainment  length,  y,  for  jets.  The  assumption  is 
made  that  for  two  different  sets  of  conditions  for  jets  (Identified  by 
subscripts  1  and  2),“ 

La/Lj  =  ya/yi  or  L2  =  ^  Lj  (l) 

where,  L2  is  the  flame  length  sought,  L,  is  the  value  obtained  from 
Fig.  5,  and  y2  and  y2  are  the  corresponding  entrainment  lengths. 
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The  calculation  of  the  entrainment  length  y  is  based  upon  the 
following  equation:18 


m  P  i 

+  !  =  0.32  f-  (-£)* 

mj  2a  PJ 


(2) 


where  m  is  the  mass  flux  of  ambient  air  which  has  been  entrained  by 
e 

the  jet  at  downstream  distance  y,  is  the  mass  flux  of  the  original 
jet  (i.e.,  before  any  mixing  has  occurred)  ,  and  p  and  p  are  the  densi- 
ties  of  the  unmixed  ambient  and  jet  streams  respectively.  In  this  report 
the  definitions  of  y  and  m  /m  are  further  restricted  to  the  special 

®  j 

values  where  M  =  m  /m  represents  a  stoichiometric  ratio  for  whatever 
®  J 

gases  are  considered.  Also  let  P  =  p  /p  .  Then  one  can  solve  equation 

e  J 

(2)  for  the  appropriate  values  of  y,  and  obtain 


Zl  = 
yi 


M,  +  1VP,; 


Substituting  values  corresponding  to  conditions  for  Figs.  4  and  5,  one 
obtains 


Zi  =  0.0513  (M2  +  1)P2"^  (3) 

y  i 

An  indication  of  the  maximum  altitude  at  which  continuum  theory 
can  be  expected  to  hold  can  be  seen  from  comparison  of  ambient  mean  free 
path,  X,  with  expanded  jet  diameter,  2a,  Table  I.  If  one  employs  the 
criterion  that  the  characteristic  dimension  should  be  at  least  100  times 
the  mean  free  path,  the  following  conclusions  can  be  drawn:  for  the  jet 
as  a  simple  blunt  body,  continuum  theory  may  hold  to  an  altitude  of 
110-120  km;  for  turbulent  mixing  assuming  a  typical  eddy  size  of  one- 
tenth  the  jet  diameter  as  the  characteristic  dimension,  continuum  theory 
may  hold  only  to  around  90  km;  for  a  shock  standoff  distance  of  0.075  times 
the  jet  diameter  as  characteristic  dimension,  the  corresponding  altitude 
is  somewhere  near  90  km.  Of  course,  the  mean  free  path  in  the  shock  layer 
itself  is  much  smaller  than  ambient  (by  a  factor  of  20-30)  because  of 
the  high  density.  As  the  size  of  the  expanded  jet  depends  upon  the  size 
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and  nature  of  the  rocket  and  its  trajectory,  these  factors  must  be 
considered  in  any  estimate  of  the  altitude.  A  quite  different  crit¬ 
erion  for  the  maximum  altitude  to  which  continuum  theory  should  apply 
is  the  relaxation  time,17  l.e.,  the  mean  time  between  collisions.  For 
large  relaxation  times,  small  disturbances  will  propagate  in  both  longi¬ 
tudinal  and  transverse  modes.  When  air  density  becomes  less  than  about 
10“ 10  atmosphere  (or  10"7  mb),  undamped  small  disturbance  propagation 
tends  to  occur.  For  the  ambient  air,  this  corresponds  to  about  105  km. 

In  the  jet,  because  the  density  is  higher  due  to  higher  pressure  and 
lower  temperature,  the  corresponding  altitude  is  higher. 

At  the  present  time  we  lack  field  data  in  "flame''  lengths  for 
various  altitudes  in  the  range  considered  here,  and  so  have  no  way  of 
comparing  with  observed  trends  of  flame  lengths  versus  altitude.  Esti¬ 
mates  can  be  made  of  available  times  for  reaction.  Dividing  the  calculated 
lengths  by  the  corresponding  rocket  velocities  provides  maximum  times  of 
mixing;  but  the  reaction  time  must  be  much  smaller  (say  not  more  than 
one-tenth  as  great)  in  order  not  to  increase  the  flame  length,  which  is 
calculated  on  the  assumption  of  a  reaction  time  small  compared  to  the 
mixing  time.  Rough  values  of  reaction  times  (dividing  the  mixing  time 
by  10)  are  then  found  to  vary  from  0.23  sec  at  80  km  to  1.5  sec  at  140  km. 
These  times  can  be  used  at  least  tentatively  in  evaluating  the  signifi¬ 
cance  of  various  electron  production  and  removal  reactions  at  different 
altitudes . 

B.  Negative  Ion  Formation 

A  previous  study*  of  ionization,  including  negative  ions  in  the 
afterburning  regions  of  rocket  exhaust,  described  qualitatively  some 
of  the  effects  which  one  would  expect  to  find  there.  However,  it  did 
not  consider  many  of  the  finer  details  of  the  ionic  processes  such  as 
the  possible  importance  of  negative  ion  formation  as  a  mode  of  electron 
removal.  Some  of  these  details  will  be  developed  in  the  following,  and 
from  the  results  we  shall  make  order  of  magnitude  estimates  of  negative 
ion  concentrations,  required  electron  production  rates  for  a  given  level 
of  ionization,  etc. 
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In  order  to  compute  the  concentrations  of  negative  ions  present  in 
the  plume  one  must  know  or  at  least  have  a  reasonable  model  of  the 
rocket  exhaust  composition.  Of  the  various  estimates  which  have  been 
made  we  shall  employ  that  proposed  by  Boynton  and  Neu13  for  a  liquid 
fueled  rocket  (Table  II).  In  addition  we  shall  assume  that  the  plume 
constituents  mix  thoroughly  with  the  ambient  atmosphere  and  that  the 
pressure  and  temperature  in  the  afterburning  region  are  equal  to  the 
ambient  pressure  and  temperature,  given  by  the  ARDC  1959  Model  Atmosphere. 
Because  of  only  partial  mixing  of  the  rocket  exhaust  with  the  atmosphere 
and  the  presence  of  shockwaves  which  result  in  higher  pressures  and 
temperatures  in  regions  of  the  plume,  these  assumptions  are  not  completely 
justified.  However,  in  view  of  the  very  crude  estimates  of  rate  constants 
which  we  shall  introduce  later,  it  is  considered  unlikely  that  errors 
associated  with  these  assumptions  will  be  greater  than  those  associated 
with  the  rate  constants.  The  lifetime  of  the  afterburning  region  is 
estimated  to  be  of  the  order  of  1  sec. 

Electron  removal  processes  are  basically  of  two  types,  recombination 
and  attachment.  Before  discussing  the  former  we  shall  consider  the 
formation  and  destruction  of  negative  ions.  The  most  important  such 
processes  are 


Radiative  Attachment 

02  +  e  -*  0a  +  hv  (al) 

0  +  e  -»  0  +  hv  (a2) 

H  +  e  ->  H  +  hv  (a3) 

Colli slonal  Attachment 

0a  +  02  +  e  -»  0a  +  0a  (bl) 

0  +  02  +  e  -»  0  +02  (b2) 
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Photodetachment 


Of  +  hv  — ♦  Of  +  e 

(cl) 

0  +  hv  -*  0  +  e 

(c2) 

H  +  hv  -*  H  +  e 

(c3) 

OH-  +  hv  -■»  OH  +  e 

(c4) 

MOf  +  hv  — *  NO j  +  e 

(c5) 

Of  +  hv  -*  Of  +  e 

(c6) 

Colli slonal  Detachment 

Of  +  Of  — *  Of  +  Of  +  e 

(dl) 

Of  +  0  >  Of  +  e 

(d2) 

0~  +  0  -»  Of  +  e 

(d3) 

H~  +  H  -*  Ha  +  e 

(d4) 

OH_  +  H  ->  H20  +  e 

(d5) 

N  +  NOf  -*  Na  +  0*  +  e 

(d6) 

0  +  H~  OH  +  e 

(d7) 

0  +  Of  — *  Of  +  Of  +  e 

(d8) 

H  +  Of  — ►  OH  +  Of  +  e 

(d9 ) 

Charge  Transfer  and  Ion-Atom 

Interchange 

Of  +  0  — *  Oj  +  0~ 

(el) 

Of  +  0  +  Of  — *  Of  +  Of 

(e2) 

Of  +  N  +  M  NOf  +  M 

(®3) 

O"  +  OH  0  +  OH" 

(e4) 

H-  +  OH  ->  H  +  OH" 

(e5) 

Oj  +  OH  Of  +  OH" 

(e6) 

Of  +  H  — »  Of  +  H 

(e7) 

0  +  Of  +  Of  -*  Of  +  Of 

(e8) 

h"  +  0  -*  o"  +  H 

(e9) 
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Reaction  (e4)  is  energetically  permitted  only  if  the  electron  affinity 
of  0"  is  less  than  1.8  ev.  Although  the  photodetachment  experiments  of 
Branscomb  and  co-workers11*  indicated  that  this  condition  is  satisfied, 
more  recent  work  by  Schulz30  (1962)  yields  an  electron  affinity  of  2  ev. 
Because  we  are  not  yet  certain  which  value  applies  to  the  0'  in  the  ion¬ 
osphere,  we  do  not  know  whether  the  process  is  possible  or  not;  for 
present  purposes  we  shall  assume  that  it  is.  Some  of  the  rate  constants 
associated  with  the  reactions  listed  above  have  been  measured.  Estimates 
of  others  are  based  on  typical  experimental  values  associated  with  similar 
reactions.  Temperature  dependence  has  been  neglected  because  of  the 
probable  absence  of  activation  energies;  the  error  so  introduced  is  ex¬ 
pected  to  be  less  than  other  uncertainties  in  the  rate  constants.  The 
experimental  results  and  estimates  are  shown  in  Table  III. 


Table  III 
RATE  CONSTANTS 


Process 

Rate 

Constant 

Value 

Reference 

al 

Pi 

10“ 18  cm*/ sec 

Branscomb  et  al.al  (1958),  using  photo- 
detachment  data  and  principle  of  detailed 
balancing 

a2 

Pa 

10”18  cm*/ sec 

Estimated  by  Poppoff  and  Whitten,** 

(1962)  using  photodetachment  data  and 
principle  of  detailed  balancing  (Burch 
et  al . , 18  1958) 

a3 

Ps 

10-16  cm*/ sec 

Same  as  a2,  using  data  of  Smith  and 

Burch**  (1959) 

bl 

ki 

10”*°  cm*/sec 

Chanin  et  al.24  (1959) 

b2 

ka 

10” 31  cme/sec 

Estimate,  based  on  idea  that  bl  probably 
involves  resonance  and  thus  should  be 
perhaps  an  order  of  magnitude  faster 

cl 

pj 

0.1  sec”1 

Estimate  (Poppoff  and  Whitten,**  (1962)) 

c2 

Pa 

1.4  sec”1 

Branscomb  et  al..2  1  (1958) 

c3 

P* 

8  sec”1 

Computed  from  results  of  Smith  and 

Burch**  (1959) 

i  c4 

P« 

10”*  sec”1 

Estimate 

c5 

Pa 

10“*  sec”1 

Estimate 

c6 

Pe 

10“*  sec”1 

Estimate 

■-  ■  ■  ..  .  . . .  ■  ■  —  ,  -  ■  1  .. 
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Table  III  (Concluded) 


Process 

Rate 

Constant 

Value 

Reference 

dl 

7i 

10~2°  cm8/ sec 

Phelps  and  Pack28  (1961) 

d2 

7% 

10" 18  cm8/ sec 

Estimate;  could  well  be  an  order  of 
magnitude  larger 

d3 

7a 

10”18  cm8/ sec 

Same  as  d2 

d4 

74 

10“18  cm8/ sec 

Same  as  d2 

dS 

7s 

10- is 

Same  as  d2 

d6 

7« 

10- is 

Same  as  d2 

d7 

7? 

10" 18 

Same  as  d2 

d8 

7s 

10- is 

Same  as  d2 

d9 

7s 

10- is 

Same  as  d2 

el 

10“ 11 

Estimate;  rate  constants  for  positive 
ion  charge  transfer  processes  in  atmos¬ 
phere  have  been  found  to  be  of  order 

10” 11  cm8/ sec,  Dickinson  and  Sayers, 26 
(1960) 

e2 

ka 

10” 88  cm®/ sec 

Estimate;  maybe  as  large  as  10”81  cm®/sec 

e3 

k3 

10”  81  cm®/ sec 

Estimate 

e4 

10”11  cm8/ sec 

Same  as  el 

e5 

10” 11  cm8/ sec 

Same  as  el 

e6 

k« 

10“11  cm8/sec 

Same  as  el 

e7 

k7 

10”11  cm8/ sec 

Same  as  el 

e8 

k8 

10”8*  cm8/sec 

Same  as  e2 

e9 

k9 

10”11  cm8/ sec 

Same  as  el 

For  a  given  electron  density  [e]  this  set  of  constants  can  be  used  in 
conjunction  with  the  corresponding  rate  equations  and  the  concentrations 
of  neutral  species  given  in  Table  IV  to  compute  the  concentrations  of 
the  various  species  of  negative  ions.  Actually  we  also  need  to  know  the 
ion-ion  recombination  coefficients  a  and  the  positive  ion  concentrations 
[X  ].  For  our  purposes,  however,  neglect  of  this  term  will  not  materially 
affect  the  final  results.  After  carrying  out  the  computations  for  the 
negative  ion  concentrations  after  1  sec,  we  obtain  the  results  presented 
in  Table  V. 
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Table  IV 


I 


CONCENTRATIONS  OF  NEUTRAL  SPECIES  (cm-3) 
Mlnzner , 13  (1959);  Handbook  of  Geophysics , 1 1  (i960); 
Boynton  and  Neu,13  (i960) 


Altitude 

[oa] 

[Kt] 

B 

[N] 

[OH] 

[H] 

80  km 

1014 

4  x  1014 

1012 

107 

1013 

10“ 

100  km 

1013 

6  x  10,a 

3  x  1012 

107 

10“ 

5  x  10® 

130  km 

1010 

1011 

10“ 

107 

10® 

5  x  107 

Table  V 

NEGATIVE  ION  CONCENTRATIONS11 


Altitude 

[o;j 

| 

[h~] 

[OH~] 

[NOj  ] 

[o;j 

.  b 
day 

80  km 

200 

350 

50 

10  s 

100 

«  1 

100  km 

0.3 

10  4 

0.1 

104 

«  1 

«  1 

130  km 

«  1 

10s 

«  1 

10 

«  1 

«  1 

Nightb 

80  km 

200 

4  x  10s 

104 

2  x  10® 

100 

«  1 

100  km 

0.3 

2  x  104 

0.1 

2  x  104 

«  1 

«  1 

130  km 

«  1 

10  3 

«  1 

10 

«  1 

«  1 

*We  are  assuming  that  [e]  =  107  cm-3,  Smith,3  (1961) 

kphotodetachment  ceases  at  sunset  and  thus  causes 
the  negative  ion  concentrations  to  rise 
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From  the  results  of  Table  V  we  can  conclude  that  [OH  ]  and  [0  ] 
are  expected  to  be  the  most  Important  negative  ions  in  the  rocket  plume. 
Although  an  increase  in  the  charge  transfer  rate  constants  will  have 
the  effect  of  increasing  the  relative  concentrations  of  N0a  the  net 
effect  will  be  small  since  it  is  extremely  unlikely  that  the  true  values 
of  the  k's  will  be  more  than  about  one  order  of  magnitude  larger  than 
those  listed  in  Table  III.  Increasing  the  values  of  the  collisional 
detachment  coefficients  y  will  have  an  important  effect,  particularly 
at  night,  by  reducing  the  total  negative  ion  concentrations,  perhaps  by 
an  order  of  magnitude  below  those  given  in  Table  V. 


Recombination  can  occur  by  either  of  two  mechanisms,  dissociative, 


e-g., 


Oj  +  e  -»  O  +  O 


(5a) 


or  ion-ion,  e.g. , 

O,  +  o;  -  Oa  +  02  ,  (5b) 

which  are  related  to  an  effective  recombination  coefficient  a  by 

eff 

the  equation 


“eff  ‘  “D  *  1  “l  <6> 

where  X  is  the  negative  ion  to  electron  concentration  ratio.  For  some 
of  the  ions  of  interest  in  the  plume  (e.g.,  CHO+,  Hs0+)  is  probably 
of  the  order  10“8  cm*/ sec  while  is  expected  to  be  considerably 
smaller,  i.e.,  ~  10”8  cm*/sec.*7  Using  the  estimates  of  negative  ion 
concentrations  in  Table  V,  we  see  that  the  dissociative  mode  is  much 
the  faster  recombination  mechanism  for  the  polyatomic  ions  and  we  can 
neglect  negative  ion  formation  as  a  deionization  process  in  the  after¬ 
burning  region  except  in  the  case  of  alkali  metal  ions  which  cannot 
undergo  dissociative  recombination.  Because  of  the  unavailability  of 
the  dissociative  mode  the  effective  recombination  coefficient  for  these 
species  is  expected  to  be  much  smaller,  probably  of  the  order  of  10”10  cm3 
sec”1 . 
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The  dominant  ionization  mechanism  in  the  plume  is  thought  to  be 
chemionization,  specifically* 

CH  +  0  -*  CHO+  +  e  (7) 

The  dissociative  recombination  coefficient  of  a  complex  ion  of  this  type 
is  probably  of  the  order  10“®  cm3/ sec,  thus  requiring  an  electron  production 
rate  of  ~  10*  cm-3  sec-1  to  produce  electron  concentrations  of  the  order 
107  cm-3.  Assuming  that  the  rate  constant  of  process  (7)  is  of  the  order 
10“ 13  cm3  sec"1,  which  is  reasonable  for  two-body  chemical  reactions,  one 
finds  that  concentrations  of  CH  of  ~  10®  cm*3  are  required.  This  is 
equivalent  to  ~  0.3  ppm  at  80  km,  ~  3  ppm  at  100  km,  and  ~  5000  ppm  at 
130  km.  The  estimates  at  80  and  100  km  appear  reasonable  but  that  at 
130  km  is  undoubtedly  much  too  large.  We  conclude  that  at  an  altitude 
in  the  region  of  about  150  km  the  electron  density  in  the  rocket  plume 
probably  approximates  the  ambient  concentration. 

If  all  of  the  positive  ions  were  alkali  metal  ions,  e.g.,  Na+,  K+, 
a  much  smaller  ionization  rate  would  be  required  to  maintain  a  given 
electron  concentration;  this  is  due  to  the  smaller  effective  recombina¬ 
tion  coefficients  previously  discussed.  Whether  or  not  the  ionization 
due  to  alkali  metals  in  liquid  fuel  rockets  is  large  enough  to  rival 
the  chemionization  process,  equation  (7),  is  not  known.  If  it  is, 
negative  ions  must  be  of  great  importance  to  the  ionization  properties  of 
the  plume.  This  point  deserves  further  investigation.  Alkali  metals 
are  present  as  Impurities  in  fuels,  and  seme  atoms  will  be  thermally 
ionized  in  the  rocket  engine.  There  is  also  the  possibility  of  photo- 
ionization  at  high  altitudes.  Also  formation  of  Cl  may  be  important 
for  aluminum  perchlorate  propellants,  especially  at  lower  altitudes 
where  Cl  concentration  is  much  greater  than  0  concentration. 

C.  Thermal  Ionization  of  Alkali  Metals 


A  high  concentration  of  electrons  may  be  present  in  the  exhaust 
gases  emerging  from  a  rocket  motor.  A  principal  source  of  these  electrons 
is  thermal  ionization  of  species  of  low  ionization  potential,  such  as 
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alkali  metal  contaminants  in  the  propellant.  High  levels  of  ionization 
presumably  are  attained  at  the  high  temperature  of  the  combustion  chamber 
some  of  this  ionization  may  be  frozen  during  nozzle  expansion  and  during 
subsequent  expansion  of  the  plume  to  ambient  pressure. 

Electrons  may  also  be  produced  by  interactions  of  the  plume  with 
the  external  environment.  If  processes  such  as  impact  ionization, 
photoionization,  and  chemionization  are  to  be  significant  sources  of 
ionization,  they  must  produce  electrons  in  greater  concentrations  than 
the  background  electron  density  due  to  thermal  ionization. 

We  shall  formulate  approximate  expressions  for  the  fractional 
ionization  of  potassium  in  the  exhaust  and  for  the  resultant  electron 
density  in  the  expanded  plume,  due  to  thermal  ionization  in  the  rocket. 

We  can  then  estimate  electron  densities  in  the  exhaust  from  this  source 
under  various  conditions.  Such  estimates  also  provide  a  criterion  for 
evaluating  the  significance  of  other  sources  of  ionization  in  the  exhaust 
plume. 


The  equilibrium  constant  K  for 

eq 

A  £  A+  +  e  (1) 

can  be  calculated  by  the  Saha  equation,  which  for  alkali  metals  is 


eq 


[A+  ]n 

TaT 


-  =  2.41  x  1015  T3/ 3  e~I/kT  (2) 


where  T  is  the  temperature  (°K),  I  is  the  ionization  potential  of  A,  k 
is  Boltzmann's  constant,  ng  is  the  electron  density  (cm-3),  and  the 
brackets  denote  concentrations  per  cm3. 
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In  the  absence  of  any  other  process  producing  or  removing  electrons, 

[A+]  =  n  ,  and  the  electron  density  is  given  by 
e 

n  =  {[A]  2.41  x  1018  T3'8  e-I/kT)*  (3) 

e 

and  the  fractional  Ionization  of  A  by 

2.41  x  10,s  T*/a  e"I/kT,J 

The  concentration  of  A  is 

[a1-'aH  <5) 

where  f.  is  the  number  fraction  of  A  in  the  gas  and  p  is  the  total 
A 

pressure , 


+] 

TaT 


We  assume  that  the  ionization  of  potassium  at  the  exit  plane  of 
the  motor  can  be  described  as  equilibrium  ionization  at  some  temperature 
T^.  This  effective  ionization  temperature  may  not  be  the  actual  tempera¬ 
ture  at  the  exit  plane,  since  ionization  may  have  been  frozen  at  some 
earlier  stage  of  nozzle  expansion.  Our  ability  to  predict  T^  is  therefore 
limited,  but  T^  will  lie  between  throat  and  exit  plane  temperatures. 

The  equilibrium  constant  for  potassium,  calculated  from  equation  (2) 
for  several  values  of  T^ ,  is  given  in  Table  VI.  We  combine  equations  (4) 
and  (5)  to  obtain  an  expression  for  the  fractional  ionization  of  potassium, 


[K+J 

Tier 


K  .  K 

-  i  -  i 

~  (TkT)  -  ( 


1  iv  k  T  i 

_  /  eq  1  \t 


f  p 

K  i 


(6) 
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Table  VI 


EQUILIBRIUM  CONSTANT  FOR 
K  <=*K+  +  e  FOR  VARIOUS 
TEMPERATURES 


T 

K 

*i 

eq 

1500 

4.0  x  10s 

2000 

2.7  x  10® 

2500 

5.7  x  1011 

3000 

2.1  x  101S 

The  fractional  ionization  thus  depends  on  the  potassium  concentra¬ 
tion  and  the  effective  ionization  pressure  as  well  as  on  Tt .  In  Table  VII 
are  given  values  of  [K+]/(Kj  calculated  for  typical  conditions,  fK  =  100 
parts  per  million  by  number,  and  Pt  *  5  psia.  For  conversion  to  other 
conditions,  note  that  [K+J/[K]  scales  as  (f^  P^) 


Table  VII 

FRACTIONAL  IONIZATION  OF  POTASSIUM  AT  VARIOUS 
IONIZATION  TEMPERATURES 


Ti 

1500 

2000 

2500 

3000 

[K+]/[K] 

5  x  10"® 

5  x  10"* 

.  .  - 

1 

o 

fH 

X 

5  x  10-1 

The  calculated  values  in  Table  VII,  as  would  be  expected,  show  a 
general  tendency  toward  increasing  fractional  ionization  as  T^  increases. 
These  values  should  not  be  compared  in  detail  because  of  the  arbitrary 
assumption  that  P^  =  5  psia.  However,  an  alternative  assumption  of 
constant  number  density  rather  than  constant  pressure  would  change  the 
relative  values  at  the  highest  and  lowest  temperatures  in  Table  VII  by 
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a  factor  of  only  v2  .  We  can  conclude  that  any  other  process  of  potassium 
ionization  which  produces  a  fractional  ionization  of  less  than  10~s  is 
probably  not  significant  in  comparison  with  thermal  ionization. 

If  we  assume  the  fractional  ionization  at  the  exit  plane,  as  estimated 
in  Table  VII,  is  frozen  during  plume  expansion  to  ambient  pressure,  then 
we  can  estimate  typical  electron  densities  in  the  plume  using  equation  (3). 
The  value  of  [K]  for  the  expanded  plume  is  given  by 

fK 

(Kj  =  ^  Pe  x  8  x  1021  .  (8) 

e 

Approximate  values  for  Pg ,  the  ambient  pressure  in  atmospheres,  are 

given  for  various  altitudes  in  Table  VIII.  These  data  are  slightly 

different  from  those  used  in  the  first  section,  but  the  differences  are 

not  significant  for  our  purposes.  The  temperature  of  the  expanded  plume, 

T  ,  can  be  calculated  for  a  particular  motor  at  a  particular  altitude 
© 

for  a  specific  expansion  model  (such  as  assumed  isentropic  expansion). 

We  shall  assume  T  =  100°K  in  order  to  calculate  typical  values  of  n  . 

0i  e 

Note  that  n  scales  as  (T  )  .  The  calculated  electron  densities  in 

e  e 

the  expanded  exhaust  at  various  ionization  temperatures  and  altitudes 
are  given  in  Table  IX;  as  before,  we  have  assumed  f  =  100  ppm  by 

IV 

number  and  P^  =  5  psia. 


Table  VIII 

AMBIENT  PRESSURE  AT  VARIOUS 
ALTITUDES 


Altitude 

(km) 

Pressure 

(atm) 

200 

10"® 

150 

3  x  10"® 

100 

10-7 

80 

1 

o 

rH 

60 

2  x  10-4 

40 

3  x  10"* 
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Table  IX 

CALCULATED  ELECTRON  DENSITIES  (cm-*)  IN  THE  EXHAUST 
AT  VARIOUS  IONIZING  TEMPERATURES  AND  ALTITUDES 


Altitude 

(km) 

Ti 

1500 

2000 

2500 

3000 

200 

2.8  x  10* 

3.6  x  104 

5.6  x  10* 

3.8  x  10® 

150 

1.1  x  10s 

1.1  x  10* 

1.7  x  10® 

1.2  x  107 

100 

3.8  x  104 

3.6  x  10® 

5.6  x  107 

3.2  x  10® 

80 

3.8  x  10* 

3.6  x  10® 

5.6  x  10® 

3.8  x  10l° 

60 

7.5  x  107 

7.2  x  10® 

1.1  x  10“ 

7.7  x  10“ 

40 

1.1  x  10® 

1.1  x  1011 

1.7  x  1012 

1.2  x  1013 

At  a  given  altitude,  ne  increases  with  ionization  temperature,  although 

the  calculated  values  are  dependent  on  the  assumed  P^  ,  as  before.  Since 

we  have  assumed  a  constant  total  number  of  electrons  in  the  plume  for  a 

particular  T^ ,  the  electron  density  naturally  decreases  with  ambient 

pressure  as  the  altitude  increases.  However,  since  T  also  decreases 

e 

with  pressure,  the  rate  of  decrease  of  ng  with  increasing  altitude  is 
not  as  large  as  indicated  by  the  values  in  Table  IX. 

Many  assumptions  and  approximations,  in  addition  to  those  already 
mentioned,  are  implicit  in  the  above  calculations.  Some  of  these  will 
now  be  specified  in  more  detail.  Sources  of  electrons  other  than 
potassium  have  been  neglected.  With  its  relatively  low  ionization 
potential,  potassium  in  the  amounts  typically  present  in  solid  propellants 
would  quite  reasonably  be  expected  to  dominate  thermal  ionization. 

Electron  sinks  other  than  K+  have  been  ignored  in  the  formulation  of 
equilibrium  ionization.  The  most  serious  effect  of  this  type  is  probably 
due  to  the  presence  of  electron-attaching  species  such  as  Cl.  The  effect 
of  equilibria  such  as  Cl  Cl  +  e  can  be  calculated  in  theory.  The 

result  would  be  lower  values  for  n@  and  higher  values  for  [K+]/[k]  than 
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given  in  tables  VI -IX.  Since  equilibrium  shifts  toward  detachment  with 
increasing  temperature,  the  assumption  of  neglecting  attachment  is  better 
at  higher  temperatures. 

It  has  been  assumed  that  the  potassium  is  available  as  atoms  in 
the  vapor  phase.  One  should  really  consider  dissociation  equilibria  for 
the  potassium  compounds  in  the  fuel  and  other  chemical  sinks  for  potassium, 
such  as  KOH  K  +  OH .  Neglecting  such  effects,  as  in  the  above  calcula¬ 
tions,  can  be  regarded  as  merely  a  limitation  on  our  ability  to  estimate 
[K]  accurately.  However,  this  should  be  given  careful  consideration  in 
future  studies. 

The  assumption  that  ionization  is  frozen  during  plume  expansion 
seems  reasonable  for  conditions  of  interest.  If  K+  and  electrons  undergo 
3-body  recombination  with  a  rate  constant  of  about  10“27  cm8  sec-1,  then 
the  time  required  to  reduce  the  electron  density  by  one-half  is  given  by 

tj  =  (10-27[M]ne)-1  (9) 

If  n@  =  1010  cm-3  and  [M]  =  1017  cm'1  (0.01  atm  at  300°K),  then 
tj  ~  1  second,  a  time  which  is  long  enough  compared  with  expansion  time 
for  recombination  to  be  neglected. 

We  concluded  that  the  amount  of  potassium  ordinarily  contained  in 
solid  propellants  may  produce  electron  densities  in  the  exhaust  plume 
at  low  altitudes  which  are  sufficiently  high  to  cause  an  attenuation 
problem.  Nevertheless  because  of  the  uncertainties  mentioned  above 
one  cannot  make  a  positive  prediction.  As  the  exhaust  temperature 
increases,  this  effect  persists  to  higher  altitudes.  The  lower  tempera¬ 
ture  in  the  combustion  chamber  and  the  lower  concentration  of  alkali 
metals  in  liquid  fuel  rockets  lead  to  considerably  lower  calculated  values 
for  thermal  ionization. 

The  purpose  of  these  approximate  calculations  is  not  only  to 
indicate  conditions  under  which  thermal  ionization  leads  to  large 
electron  densities  ,  but  also  to  provide  a  criterion  for  estimating  the 
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importance  of  other  ionization  mechanisms.  Obviously,  only  those  reactions 
which  can  produce  electrons  in  larger  concentrations  than  the  background 
thermal  ionization  can  have  an  ^important  effect  on  the  ionization  level 
of  the  exhaust. 

For  example,  consider  a  general  process: 

X*  +  K  -»  X  +  K+  +  e  (10) 

where  X*  is  an  excited  species  whose  excitation  energy  is  greater  than 
the  ionization  potential  of  potassium.  If  this  process  is  to  have  a 
significant  effect  on  the  electron  density,  the  concentration  of  X*  must 
obviously  be  larger  than  the  background  electron  density,  or  its  rate 
of  production  in  the  exhaust  must  be  much  larger  than  the  required 
electron  production  rate.  The  electron  production  rate  by  this  process 
must  exceed  the  recombination  rate  if  an  increase  in  electron  density 
is  to  occur.  As  an  upper  limit  to  the  recombination  rate  assume  dissoci¬ 
ative  recombination  with  a  recombination  coefficient  of  10-®  cm3  sec-1 
can  occur.  Then,  if  the  rate  constant  for  reaction  (10)  is  about  10-12  cm3 
sec-1,  we  require  that 

10-12[X*][JC]  >  10-8[K+]2  .  (11) 

If  [K+]/[K]  ~  10-4  and  [K ]  -  1010  cm-3,  then  equation  (11)  indicates 
that  a  concentration  [X*]  >  10®  cm-3  is  sufficient  to  increase  the 
electron  density. 

A  more  convenient  method  for  applying  this  criterion  may  be  to 
define  a  characteristic  time  for  electron  loss  (by  recombination  or 
diffusion)  in  the  plume.  Then,  if  the  electron  production  rate  by  a 
particular  mechanism  can  be  estimated,  the  total  ionization  per  unit 
volume  during  this  characteristic  time  can  be  compared  directly  with 
the  background  thermal  ionization.  For  example,  the  rate  of  photoioniza¬ 
tion  of  potassium  by  solar  radiation  is  — j—  sec-1,  where  F  is  the 
total  flux  per  cm2  per  sec  in  the  appropriate  wavelength  region,  Q  is 
an  averaged  photoionization  cross  section,  and  I  is  the  ionization 
potential.  If  the  background  fractional  ionization  of  potassium  is 
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[K+]/[K],  the  time  required  for  photoionization  to  increase  the  ioniza- 
tion  significantly  is  .  This  must  be  less  than  the  characteristic 

time,  if  photoionization  is  to  be  an  important  process. 

D.  Photoionization 

The  importance  of  photoionization  effects  in  missile  exhaust  plumes 
is  being  explored.  A  preliminary  look  at  this  problem  indicates  that 
sufficient  solar  energy  in  the  ultraviolet  wavelengths  may  penetrate  to 
altitudes  at  which  rocket  motors  are  still  burning  or  outgasslng.  The 
study  is,  of  course,  restricted  to  those  photon  energies  that  are  capable 
of  photoionization. 

Table  X  summarizes  estimates,  based  on  rocket  data, 28 > 29 > 50  of 
altitudes  at  which  1/e  of  the  solar  flux  (of  given  wavelength)  incident 
to  the  top  of  the  atmosphere  will  penetrate.  The  maximum  ionization 
potentials  that  can  be  supplied  at  the  1/e  penetration  altitude  by  the 
photons  are  also  listed,  along  with  estimated  fluxes  and  an  indication 
of  some  common  substances  that  could  be  ionized  at  or  above  the  indicated 
altitudes.  Except  for  Lyman-a  of  hydrogen,  the  values  are  averages  and 
do  not  represent  specific  strong  line  emission. 

Calculation  of  ionization  rates  is  rather  lengthy  and  requires 
compilation  of  considerable  data  because  the  spectral  distribution  of 
radiation  intensity  varies  with  altitude,  and  both  the  spectral  radiation 
intensity  and  the  photoionization  cross  sections  vary  (often  in  a  rather 
complicated  manner)  with  wavelength.  However  it  may  be  of  interest  to 
consider  some  examples  of  rate  of  ionization  per  unit  wavelength  interval 
for  NO,  K,  and  CO,  using  corresponding  absorption  cross  sections  of 
2  x  10"*8  cm2,  6  x  10”21  cm2,  and  1.6  x  10"17  cm2,  respectively . 31 ' 22 
These  cross  sections  correspond  to  the  wavelengths  attenuated  to  1/e 
of  the  initial  flux. 
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Table  X 

ESTIMATES  OF  PHOTOIONIZATION  PARAMETERS 


Solar  Radiation 
Attenuated  to 
l/e  at  Altitude 
(km) 

Wavelength 

of 

Attenuated 

Radiation 

(A) 

Photon 
Energy  of 
Radiation 

(ev) 

Examples  of 
Substances  Ionized 
at  or  above  l/e 
Penetration  Altitude 

Approximate  Average* 
Spectral  Radiation 
Intensi ty  l/e 
Penetration  Altitude 
(ergs/cm2  sec  A) 

23 

3000 

Cs 

16.5 

30 

2800 

K 

8.5 

40 

2600 

4.8 

5.2 

40 

2400 

5.2 

Na 

2.5 

33 

2200 

5.6 

Li 

2.3 

42 

2000 

6.2 

Al,  Ca 

0.51 

82 

1800 

6.9 

0.14 

110 

1600 

7.8 

Ni,  Mg 

0.023 

112 

1400 

8.9 

Fe,  Co 

1.9  x  10" 3 

=»75 

La 

NO 

1.87 

70-110 

1200 

10.3 

4.2  x  10“2 

95 

1000 

12.4 

C21  ,  NOa  ,  NHj, 

Paraf i nes 

1.3  x  10*3 

130 

800 

15.5 

H,  0,  Cl,  N,  H„,  Na, 
Cl3,  CO,  OH,  HaO,  CO 

ch4 

9  x  10-4 

1 

162 

600 

20.6 

F 

7.3  x  10-4 

150 

400 

31 

10*  3 

110 

44-60 

300-205 

10~2 

100-115 

8-20 

1500-620 

2  x  10“3 

65-100 

2-8 

6200-1500 

10*  3 

*Not  uniformly  averaged  some  50A,  16A,  12A  increments. 
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Examples  of  ionization  rates  that  may  be  possible  to  obtain  at 
the  l/e  penetration  altitudes  are  as  follows. 


1.  NO:  Can  be  ionized  by  radiation  shorter  than  approximately 
1200A  (Lyman  Q  would  be  the  major  radiation  for  this  reaction). 

Two  ion  pairs  cm-3  sec-1  A-1  would  be  produced  per  107  cm-*  NO 
molecules . 

2.  K:  Four  ion  pairs  cm“3  sec-1  A-1  would  be  produced  per  10®  cm-3  K 
atoms . 

3.  CO:  Five  ion  pairs  cm-3  sec-1  A-1  would  be  produced  per 
1010  cm-3  CO  molecules. 


During  the  next  quarter,  attempts  will  be  made  to  determine  the  order 
of  magnitude  of  concentrations  of  ionizable  constituents  that  can  be 
expected  to  exist  in  missile  plumes  at  various  altitudes.  If  consider¬ 
ation  of  these  data  still  indicates  that  photoionization  may  be  important, 
a  more  thorough  study  of  the  problem  will  be  planned.  In  the  meanwhile, 
compilation  of  radiation  flux  and  photoionization  cross  section  data 
will  continue. 
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